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Abstract 

We present the general lattice model for a multi-component atomic Bose- 
Einstein system in an optical lattice. Using the model, we analytically study 
the quantum phase transition between Mott insulator and superfluid. A 
mean-field theory is developed from the Mott insulator ground state. When 
the inter-species interactions are strong enough, the Mott insulator demon- 
strates the phase separation behavior. For weak inter-species interactions, the 
multi species system is miscible. Finally, the phase diagram is discussed with 
the emphasis on the role of inter-species interactions. The tips of the Mott 
insulator lobes do not depend on the inter-species interactions, but the latter 
indeed modify the range of lobes. 
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I. INTRODUCTION 



The study of quantum phase transitions (QPT) has attracted much interest in recent 
yearsSi. The term "quantum" is used to emphasize that it is quantum fluctuations that play 
a vital role in driving the transition from one phase to another. In contrast, the usual ther- 
modynamic phase transition at finite temperature is driven by thermal fluctuations which 
are experimentally controlled by tuning the temperature of the system. As temperature is 
lowered, the thermal fluctuations are suppressed and finally they are not strong enough to 
drive a finite temperature phase transition. However, this by no means implies that there 
would be no phase transition at very low temperature, since quantum fluctuations still exist 
and they may be sufficiently strong to drive a phase transition even at zero temperature. 
We call such a zero temperature phase transition a QPT, and it is experimentally accessible 
by tuning parameters of the system other than temperature. 

Several prominent examples have been extensively studied to demonstrate QPT. One 
example is quantum Hall (QH) systems, where different QH phases can be achieved by 
tuning either the magnetic field or carrier concentration^. The second examle is a network 
of Josephson junctionsi. A Josephson junction is a tunnel junction connecting two super- 
conducting metallic grains. A Cooper pair of electrons are able to tunnel back and forth 
between the grains. If the Cooper pairs can move freely from grain to grain in the network, 
the system is superconducting. However, since the grains are very small, it costs a charg- 
ing energy to move a Cooper pair to neighboring grains. When the charging energy is big 
enough, the Cooper pairs fail to propagate among the grains and the network will be in an 
insulating phase. 

A third system that exhibits QPT involves the superfluid A He. When the superfluid 
A He is absorbed in the porous media or on different substrates, the bosonic atoms in A He 
experience external forces from the other medium. When the interactions between atoms 
are much weaker than above external forces, the system is expected to be a superfluid. In 
the opposite limit, the superfluid phase can not be maintained, and the system will exhibit 
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a Mott insulator behavior. Thus a superfluid-Mott-insulator phase transition is expected 
to happen, if one can tune the strength of atomic interactions. Detailed discussions can be 
found in Ref. ^ by Fisher et al. The starting point is the following boson Hubbard model: 



Here at and a] correspond to the bosonic annihilation and creation operators on the i-th 
lattice site, rij = ajaj the atomic number operator on the z-th site, and £j the energy offset of 
the atom on the i-th site due to external harmonic confinement. The last term corresponds 
to the on-site repulsion between atoms, while the first term describes the tunneling of atoms 
between neighboring sites. At mean field level, starting with a strong coupling expansion, 
namely treating the hopping term as a perturbation, the system is found to have a QPT at 
the following critical value@"0 for the ratio U/J: 



where z = 2d for a d dimensional simple lattice and no is the inverse fraction of condensed 
atoms in a canonical ensemble. For instance, no ~ 5.83 for the three dimensional case. 

Experimentally, such critical point of QPT is very hard to access. Temperature is an 
annoying factor for a convincing demonstration of the QPT: The intrusion of thermal fluc- 
tuations often washes out the effects of quantum fluctuations. This makes the temperature 
window to observe the QPT small. Moreover, to make the system go cross the quantum 
critical point, we need to tune the controlling parameter carefully. In most of the studied 
cases, this is hard to manipulate. Even one can tune the parameter, the range of tunability 
is normally very small. Until very recently, in most cases only the magnetic fieldSi is the 
tunable parameter. Finally, the presence of disorder makes the observation of QPT even 
more difficult. 

Recently Ref. |9] reported the success in realizing a superfluid-Mott-insulator phase tran- 
sition in a gas of ultra-cold atoms in an optical lattice. This is a revolutionary breakthrough 
for experimental observation of a controllable QPT. They cooled the atomic gas of 87 Rb 
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down to 10 nK to realize the atomic Bose-Einstein condensation (BEC). Moreover, the BEC 
is loaded into a perfect, simple cubic, optical lattice formed by six criss-cross laser beams. By 
controlling the intensity of the laser beams, they can efficiently control the potential height 
of the above simple cubic lattice in a very large range. In addition, such a unique invention 
of the artificial lattice has the great advantage that the system is basically defect-free. By 
using this set-up, they successfully and repeatedly observed the QPT at the critical value 
given by Eq. @. Thus an ideal playground for QPT has been created in the atomic BEC 
system, which provides us an opportunity to test many theoretical predictions. 

A significant difference between the atomic BEC superfluid and the 4 if e superfluid is that 
the former allows atoms to condense with different internal states due to hyperfine splitting. 
This allows the order parameter of superfluid to possess a larger symmetry than the familiar 
£7(1)00. It is dubbed in the BEC community as spinor BEC. As pointed out by and 
many othersffl, the spinor BEC possesses a whole host of quantum phenomena that are 
absent in the scalar cases: For instance, vector and quadrupolar spin wave modes, Skymions 
and other quantum orders etc. Experimentally, one can condense different matter species 
into one single internal state and study the effects of cross-species interactions. Throughout 
this paper, we would like to call such systems as multi-component BEC systems. 

In the Mott insulator to superfluid quantum phase transition, quantum fluctuations 
and atomic interactions play a vital role. Without interactions, one has only the so-called 
band insulator. In the atomic gas, due to laser cooling technology, BEC can be realized 
simultaneously in several internal hyperfine levels^!. This makes the experimental study 
of the multi-component BEC possible. Among all the interesting physics discovered in the 
multi-component BEC, the inter-species repulsive interactions play a very important role. 
Therefore, it would be very interesting to study how inter-species interactions affect the 
transition from the Mott insulator to superfluid. 

Motivated both by experimental progress and by theoretical curiosity, we shall study 
in the present paper the superfluid-Mott-insulator transition in a multi-component BEC 
system in the presence of a periodic potential created by criss-cross laser beams. The layout 
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of the paper is the following: The boson-Hubbard model for the multi species is derived 
for the general case and some special cases as well in Section II. In the section III, we 
study the ground state and its stability in the strong coupling limit. The phase boundary 
between superfluid and Mott insulator is determined for the two-component case in Section 
IV. Finally, we summarize our results in the Section V. 



II. THE MODEL 

A. The general boson Hubbard model for a multi- component BEC 

After including the optical lattice potential, the most general model Hamiltonian for a 
multi-component boson gas can be written, in the second-quantization notations, as 

h 2 w 2 



H= d 6 x 



(ft , V \ 
~ 1 2m~ Sij + Uij ^ + V i( x )Sij)4>j( x ) ( 3 ) 

where m a is the mass of an individual atom, the indices i,j,k,l label the components of 
the atoms and the summation is assumed for repeated indices. Generically we allow the 
external potential C/jj to have a non-diagonal part in the hyperfine spin basis, in which it 
represents a Josephson-type coupling between spin components^. V{(x) denotes the optical 
lattice potential seen by atoms of species i. For the experimental configuration in Ref. |^, 
this lattice is modeled by 

V(x, y, z) = Vo(sin 2 kx + sin 2 ky + sin 2 kz), (4) 

with k the wave vector of the laser light and Vq the depth of the potential well. In the 
multi-component case, the depth Vq^ may depend on the species index i. The inter-atomic 
interactions in Eq. (|3]) have been approximated as a contact interaction in which the coeffi- 
cients gij t ki describe the strength of various elastic and inelastic collisions. 

For a single atom in the trap and the periodic potential, the energy eigenstates are Bloch 
states. In the tight-binding (TBA) limit, we can superpose the Bloch states to get a set of 



Wannier functions, which are localized on individual lattice sites. Within the single band 
approximation, we can expand the field operators in the Wannier basis as 

^i{ x ) = bniWi(x - X n ), (5) 
n 

where Wi(x — x n ) is the Wannier function around lattice site n. Using Eq. (||), the general 
Hamiltonian (|3D is reduced to a generalized boson Hubbard Hamiltonian for the multi- 
component BEC: 

(6) 

<m,n> m " m 

Here J^ n is the hopping matrix element between two adjacent lattice sites m and n. It is 
defined by 

r h 2 v 2 l 

J mn = ~ J d 3 XW*(x - X m ) [—^~ &ij + V i$ij + U ij ~ ^{ U H + U jj)) W j( X ~ X n)- (7) 

Within the Hubbard approximation, the hopping integral is lattice site independent, i.e., 
J m,n ~ J 1 "* • £ m describes the energy offset on each site due to the trap confinement. It is 
defined as 

d 3 xUii(x)\wi(x - x n )\ 2 . (8) 

(Here we only consider the lowest band in the optical lattice, whose bottom is taken to be 
the zero point for energy.) Finally, to get the on-site interactions in Eq. ([]), the Hubbard 
approximation has been used to approximate the multi-center integral as a single-center one; 
namely, we have 



Uij,ki = gij,ki J d xw*(x - x n i)w*(x - x n2 )w k (x - x n3 )wi(x - x nA ) (9) 
^ 9ij,ki / d 3 xw*(x)w*(x)w k (x)wi(x). 



The general form of the boson Hubbard model (§) contains a large number of parameters. 
In the following, we would like to discuss several special cases which might be relevant to 
experiments. The first simple case is, of course, given by Eq. ([[]) for a single component 



BEC. It is first derived and studied in the context of atomic BEC in Ref. 16. 



B. Two- component boson Hubbard model 



The second example we will discuss is the two-component BEC. Experimentally, the 
simultaneous condensation of 87 Rb atoms in the two internal states (F = 2, M — 2) and 
(F = 2, M = — 1) has been accomplished by Myatt et al.0. For this case, we discuss Bose 
condensed atoms with two internal hyperfme levels \A > and \B >. The atoms interact only 
through the following three channels: A A, BB, and AB type elastic collisions. Then our 
Hamiltonian (IBI) is reduced to 



H2= - E 



<m,n> 



(10) 



J A bl A b nA + J B b ] mB b nB + J AB bl A b nB + B.c. 

+ e mAn m A + £mBn m B J 

+ \ E U A n mA (n mA - 1) + U B n mB (n mB - 1) + U AB n mA n mB 

A similar energy-level and interaction pattern has been discussed by another group in a 
different context^. To be concrete, we focus on the situation in which the trap potential is 
diagonal in internal space, namely, £7^ only have diagonal components Uu. In this case, it 
follows from Eq. (0) that 

4 B = o. (ii) 

To get more insight into the parameters in Eq. (pi]), we have to use the explicit form of 
the Wannier functions. To do so, we notice that the optical lattice potential is sinusoidal. 
The Wannier function could be constructed as the localized one determined by the following 
eigenvalue problem: 

[^— + V(x, y, z)]<p(x, y, z) = E<j>(x, y, z). (12) 

The lattice sites are given by minima of the lattice potential V(x,y, z); around them the 
potential V is approximately quadratic: 

V(x,y,z) w ^m a u(x 2 + y 2 + z 2 ), (13) 



where u is given by 

uj = 2m a k 2 V . (14) 

(Here for simplicity, we assume that the depth of the optical potential is the same for different 
species; it is straightforward to generalize our results to the case with Vo — > Vo,i dependent on 
the species index i.) Therefore, within the single band approximation, the Wannier function 
is approximately given by the ground state wavefunction of a three-dimensional harmonic 
oscillator. Namely 



a 



iiW/u ( J-^= ) exp[-^a 2 (x 2 + y 2 + z 2 )}, (15) 



where a = y ^7^- Noting that the Wannier functions are independent of species indices 
within our approximations. Thus we can take 

Ja = Jb = J- (16) 
Finally, the on-site energy is the original inter-atomic interaction with an extra numerical 



factor J[w(x,y, z)} . Therefore, the Hamiltonian (10) can be cast into 



H 2 = - J ( b lnA b nA + ^ mB b nB + H.C. ) + ^ ( e mA n mA + E mB^mB (17) 



<m,n> 

+ ^ U An m A(n mA - 1) + U B n mB (n rnB - 1) + U AB n mA n mB 

m - 

Our Hamiltonian (|i~7D is different from the two-species boson Hubbard model proposed 
in Ref. [16], where the authors assumed that two species A and B are placed in two different 
optical lattices with a relative half-period shift. Also a drive laser has been applied to 
induce the transition between species A and B. In this situation, J A and J B in fllOD should 
be neglected since they represent the next nearest neighbor hopping. Moreover, the on-site 
mutual interactions between two species are of higher orders compared with the on-site 
interactions for the same species. In this way, we recover their Hamiltonianlil 

H 2 = - J E {bl A b nB + H.c.)+J2[ EmAnmA + E mB n mB ) (18) 



<m,n> 

U A n mA (n mA — 1) + U B n mB (n mB — 1) + U AB n mA n mB 
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C. Boson Hubbard model for spinor BEC 



Another well studied example of the multi-component BEC is the so-called spinor 
BECffl. For a system of spin / = 1 bosons, such as 23 Na, 39 K, and 87 Rb atoms, the 
form of the inter-atomic interactions is largely constrained by symmetries. In this case, the 
number of the interaction parameters are reduced to two. The interaction potential can be 
written as 

V int (x 1 , x 2 ) = (g + g^Fi ■ F 2 )S(x 1 - x 2 ), (19) 

where the parameters go and g 2 are defined by the scattering length a 2 and cto as 

Anh 2 2a 2 + a 

go = 7. — ; 20 

m a 3 
4irh 2 a 2 -a 

92 = — ■ 21 

m a 3 

To facilitate the discussion, we choose a basis to make the trapping potential in Eq. @ 
diagonal. Thus the hopping integral is non-vanishing only between the same species. In 
addition, the on-site interactions are reduced to the following eight terms: 

Hm = \ U A n mA (n mA - 1) + I Y U AB n mA n mB + ^ XX fo mo & mi&mi + H - c )- (22) 

1 mA 1 m,A^B 1 rn 

Here the species index A = 1, 0, 1(= —1). U A , U A b, and U are determined by the parameters 
go, #2 arid the Wannier functions. In particular, the parameter Uo is proportional to g 2 . 
Moreover, g 2 is determined by the difference between scattering lengths as shown in Eq. 
(^T|). For the sodium case, the difference between two scattering lengths is very small 
(0.29nm) compared with 2a 2 + a = 7.96nm. Therefore, g 2 -C go so that we can neglect the 
spin relaxation channel in the interaction terms. Namely, we set Uo = as the zeroth order 
approximation. Within this approximation, we get the following boson Hubbard model for 
a spin-1 spinor BEC, when it is loaded into the optical lattice potential: 

H 3 = - J2 ( jAb m AbnA + H.c.) + J2 tmAn mA (23) 

<m,n>A mA 

+ ^ U A n mA (n mA - 1) + - Y U AB n mA n mB - 

1 mA 1 m.A+B 



In the next section, we will discuss the possible mean field phase diagram for the superfluid- 
Mott insulator phase transition by starting with the two-component boson Hubbard Hamil- 
tonian flTop. 



III. MOTT GROUND STATE AND ITS STABILITY 

We are going to employ the strong coupling expansion to develop a mean field theory. In 
the strong coupling limit, the hopping term can be treated as a perturbation. In the zeroth 
order approximation, we ignore it for a while. The Hamiltonian is then decoupled for the 
site index. The ground state is given by the occupation number state |n^, ns >, with the 
wavef unction 



Gnd > MF ~ WmA) nA {bl B T*\U > . (24) 

m 



To get the ground state energy, we need to minimize the energy at each site (for this purpose, 
we neglect the site index in the following discussions). Namely, we need to minimize the 
energy function E(riA,nB) given by 

E(n A ,n B ) = e A n A + e B n B + ^[U A n A (n A - 1) + U B n B (n B - 1) + U AB n A n B }. (25) 

If we skip over the fact for the moment that the occupation number n A and n B are integers, 
then the condition to minimize the above energy function are given by 

U A n A + U AB n B = ~y - e A , (26) 
U AB n A + U B n B = -^--e B . (27) 

Sovling the two coupled linear equations, we get 

U B (U A - U AB ) + 2(c b U A b - e A U B ) /00 , 

UA ~ 2(U A U B - U\ B ) ' (28) 

U A {U B ~ U AB ) + 2(e A U A B - e B U A ) 

UB = 2(UaU b -UI b ) • (29) 
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Now we take care of the fact that the occupation numbers must be integer. So the actual 
numbers to minimize the energy are the two integers closest to the above n A / B . To do so, 
we can write n A / B in terms of the closest integer numbers n A / B and the decimal parts, i.e. 

n A = n° A + a, n B = n° B + (3, (30) 

where the numbers a and (5 satisfy 

- \ < « = n A - n° A < ^, (31) 
-\<P = n B -n B < 1 -. (32) 

Namely, when the parameters of the system satisfy the following conditions 

n i , Uab{Uab-Ub + 2s b )-2s a U b 
HA ~ l< 2{U A U B - U\ B ) <n ^' (33) 

o i . Uab(Uab -Ua + 2e A ) - 2e B U A 
nB_1< 2(^t/ fl - U\ B ) <ns ' (34) 

the occupation numbers {n A ,n B ) minimize the energy E(n A ,n B ). 

One loose end in above discussions is that we have assumed the minimal occupation 
numbers are non-zero. If one of the occupation numbers is zero, then it means that our 
ground state is not stable due to the mutual interactions between different species. To get 
the stability condition for the uniform ground state, we need to diagonalize the interaction 
terms, namely the following quadratic form: 

U (n A , n B ) = -{U A n 2 A + U B n 2 B + U AB n A n B )- (35) 

The eigenvalues of this quadratic form are 

n± = - [{U A + U B ) ± yf(U A - U B ) 2 + U\ B ] (36) 
1 



= - A [(U A + U B ) ± x/(U A + U B y + (U\ B - 4U A U B )\. 

Therefore, n_ may become negative; if so, the interaction manifold is saddle-like and one 
cannot really minimize the ground state energy with two non-zero occupation numbers. 
Thus, in one spatial region, one of the species must have zero occupation. In other words, 
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the ground state of the system must be phase separated, when the following condition is 
satisfied: 

U a Ub > \u AB . (37) 

This condition ( |3"TD for phase separation is analogous to that of an ordinary two-component 
BEC (without being loaded into an optical lattice )B In the case when the Wannier functions 
are the same for both species, this condition is reduced precisely to the one in the absence 
of the optical lattice. 



IV. PHASE TRANSITION TO SUPERFLUID 

In this section, we are going to present a mean-field theory based on the ground state 
developed in the preceding section. The hopping processes correspond to moving bosons 
from one site to another. This process allows bosons at different sites communicate with 
each other and finally they conspire to establish macroscopic coherence under appropriate 
conditions. In this way the system can enter a superfluid state with indefinite filling of 
bosons at each site. 

The consistent mean-field theory we shall use corresponds to the following decomposition 
of the hopping terms: 

tf m b n « <b] n >b n + b\ n <b n > - <b] n >< b n >, (38) 
= <P(bl + b n ) - 4>\ 

where 4> =<b\ n >=< b n > is the superfluid order parameter. In the case at hand, we have 
taken the order parameter to be real. In this decomposition, the higher order fluctuations 
(b^ ~ 0)(&n — 0) have been neglected. It reflects the fact that in the ground state energy 
corrections we neglect the correlation energy. Generally speaking, this process will increase 
the energy of the system; however, when the system parameters satisfy certain conditions, 
this process will not cost any energy or even will lower the energy of the system. This signals 
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the occurrence of a phase transition. Therefore, the vanishing energy correction due to the 
hopping process should give us the phase boundary. In the following, we shall determine 
the phase boundary using second order perturbation theory. 

The resulting mean-field version of the hopping Hamiltonian can be written as 



eff 



(39) 



-ZJ M^mA + Ka) + <pB{b ] mB + b mb ) 



?A + <&) 



Here z is the number of nearest-neighbor sites. Since it is a single sum over all lattice sites, 
we drop the site index from now on. 

The first order correction to the energy vanishes, due to the fact that the ground state 
is a product of number eigenstates at each site, and thus the average of an annihilation or 
creation operator is just zero. The second order correction to the energy is given by the 
following well-known expression: 

1 < g\H$*\n > | 2 



4 2) = £- 



(40) 



where |n >= \n A , n B > denotes the unperturbed state with n A and n B atoms for each 
species, respectively. Correspondingly, \g >= \n° A , n° B > is the ground state and the occu- 
pation numbers are given by Eq. (0). A straightforward calculation gives the second order 
correction to the ground state energy as follows: 

n° A + l 



t2 2 ;2 
J Z <j> A 



n A 



s A + U A (n° A - 1) + U AB n%/2 -e A - U A n° A - U AB n%/2 



(41) 



+ JV0| 



n 



B 



£b + U B (n% 
+ Jz{<f> 2 A + <f> 2 B ). 



+ 



n 



B 



1) + U AB n° A /2 -s B - U B n% - U AB n° A /2 



Therefore, the phase boundaries between the Mott insulator and the superfluid for species 
A and B, respectively, are given by the following conditions: 



1 + Jz 
1 + Jz 



+ 



n° A + l 



e A + U A (n° A - 1) + U AB n%/2 -e A - U A n° A - U AB n° B /2 



71, 



n° B + l 



e B + U B (n° B -l) + U AB n° A /2 -e B - U B n° B - U AB n° A /2 
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0. 
0. 



(42) 



Solving above equations yields 
1 



4-1 



U AB n B + U A (2n° A -1)-Jz± JU\ - 2U A Jz(2n° A + 1) + (Jz) 2 



(43) 



UabA + U B (2n° B -1)-Jz± sjUl - 2U B Jz{2n% + 1) + (Jz) 2 . 
When the chemical potential e A is in the following region 

e A < e A < 4, (44) 

the energy correction due to the tunneling events of species A is positive and thus the 
Mott-insulator is the stable ground state. A similar result is valid for species B. When the 
inter-species interactions U AB and the intra-species interactions U A / B satisfy the condition 
( |3?D , we see that the mutual interactions and non-zero occupation for the other species 
indeed modify the phase boundary for both species. When the mutual-species repulsion 
is strong enough to make the system phase separated, then the system can be viewed as 
two totally independent single species in the optical lattice. However, the most interesting 
observation from Eq. ( [43"D should be that the tips of the lobes given by condition ([H]) is 
independent of the inter-species interaction U AB and thus it is independent of the occupation 
number of the other species too! The tips of the lobes for species A and B are given by 

= 2n% B + 1 + J(2n A/B + If - 1, (45) 



\~Jz~ J c ~ A/B V ^ ' A/B 

which is the same as for the single species boson Hubbard model. The tip of the first lobe 
is given by the critical value U A / B /(Jz) = 5.83 for both species. It is worth noting that the 
above result is valid even if the hopping integrals are not equal: J A ^ Jb] the only thing we 
need to do is to scale U A / B by the corresponding hopping parameter J A / B respectively. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we have analytically studied the quantum phase transition between su- 
perfluid and Mott insulator for a multi-component BEC system in an optical lattice. The- 
oretically, this is a generalization of the well-studied case of the boson Hubbard model in 
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the condensed matter literature. Experimentally, the first beautiful observation of such a 
quantum phase transition is accomplished by loading a one-component atomic Bose-Einstein 
condensate into an artificial optical lattice and thus it would be nice to study the role of the 
inter-species interactions in the QPT for the multi-component cases. 

In the first part of the paper, we have generalized the single species boson Hubbard 
model to the multi-component case with most general interactions. To be concrete, we have 
reduced our general boson Hubbard model to the two- and three-component cases under 
appropriate conditions. 

Starting with the two-component boson Hubbard model, we developed a mean-field 
theory to study the quantum phase transition. Depending on inter-species interactions, the 
system may be in different ground states. If the repulsion between two species is not very 
strong, the two species can co-exist; namely, the system is miscible. However, if the this 
repulsion is sufficiently strong, the two species may become immiscible, and the ground state 
will demonstrate the behavior of phase separation. Namely two species of Mott insulator 
will stay in separate spatial regions. After turning on the tunneling terms, the ground state 
energy will get corrections from tunneling. We calculated the energy corrections up to second 
order and determined the boundary between the gain and loss in energy. We found that 
the inter-species interactions indeed can change the range of the parameters for the Mott 
insulator. However, the inter-species interactions can not change the position of the tips for 
the Mott insulator lobes. The phase diagram of the two-component boson Hubbard model 
also demonstrates a richer structure. From our analytical treatment, we conclude that the 
following three different phases are possible: (1) Both species A and B are in the superfluid 
phases; (2) one of the spieces is still in the superfluid phase, while the other is in the Mott 
insulator phase; (3) both species are in the Mott insulator phase. 

Finally, several further remarks are in order. (1) Above results for a two-component 
system can be directly generalized to the spinor BEC and other systems with more com- 
ponents. (2) One can check without difficulty that there is indeed an energy gap in the 
excitation spectrum of the Mott insulator. The gap at zero momentum is determined by 
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the on-site Coulomb energy Ua/b and is independent of the inter-species interaction Ua/b- 
Right at the phase transition, due to the gain of tunneling, the energy gap closes. Therefore, 
the system becomes compressible (or gapless) and therefore in a superffuid phase. (3) In 
this paper, our studies of the phase diagram for the multi species boson-Hubbard model 
have been restricted in the case where the Josephson-type tunneling term can be neglected. 
Under certain experimental conditions, such terms would be dominant and the physics is 
significantly changed. The results for this case will be published somewhere elseil. 
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FIGURES 

FIG. 1. The schematic phase diagram of the two-component boson-Hubbard model. The phase 
diagram is divided into three regions for fixed occupation numbers n^m- One region is superfiuid 
phase for both components (S.F. phase); the second region is a mixture of superfiuid phase for 
one component with Mott insulator phase for the other (S.M. phase); The third region is Mott 
insulator for both components (M.I. phase). 
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